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In this work we perform a Bayesian inference of the crust-core transition density p+ of neutron stars 
based upon the neutron-star radius and the neutron-skin thickness data within a thermodynamical 
method. The uniform and Gaussian distributions for the p; prior are adopted in the Bayesian 
approach. It has a larger probability to have values higher than 0.1 fm~? for ps as the uniform prior 
and the neutron-star radius data are used. We found that this is controlled by the curvature Ksym of 
the nuclear symmetry energy. This phenomenon will not happen if Ksym is not extremely negative, 
i.e., Ksym > -200 MeV. The obtained p+ is 0.075+9:095 fm~* at 68% confidence level when both the 
neutron-star radius and the neutron-skin thickness data are taken into account. The strongly anti- 
correlations between p; and the slope L, curvature of the nuclear symmetry energy are observed. 
The dependence of the three L-Ksym correlations predicted in the literature on the crust-core density 
and pressure is quantitatively investigated. The most probable value of 0.08 fm~* for p+ is obtained 
from the L-Ksym relation raised by Holt et al. and the larger values are preferred by the other two 


relations. 
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I. INTRODUCTION 


The determination of the crust-core transition densi- 
ty pz: in neutron stars (NSs) is important for not only 
the predictions of the bulk NS properties[1] but also the 
finite nuclei properties[2—4]. However, due to the intri- 
cate structure of the inner crust in NSs, constraining the 
transition density is still challenging. Various theoretical 
models have been used to estimate the transition densi- 
ty in the past years including the dynamical method[5— 
9], the Thomas-Fermi approximation method[10], the 
random phase approximation|4], the thermodynami- 
cal method[11-13], the Vlasov method[14], the com- 
pressible liquid-drop model{15] and the meta-modeling 
approach[16]. These methods give different prediction- 
s, i.e., p,=0.07140.011 fm? estimated in the meta- 
modeling approach[16], ;=0.0955+0.0007 fm? ob- 
tained by comparing with the excitation energies of gi- 
ant resonances, energy-weighted pygmy dipole strength, 
and dipole polarizability data using the relativistic nu- 
clear energy density functionals[2], p;=0.04~0.065 fm~* 
limited by using the EOS including the momentum- 
dependent interaction of neutron-rich nuclear matter 
constrained by the isospin diffusion data in heavy-ion 
reactions|8], p:=0.069~0.098 fm~? from the thermody- 
namical method[13] and p,=0.058~0.092 fm~* with the 
Thomas-Fermi method[10] and so on. These predicted by 
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the dynamical method are usually smaller than that from 
the thermodynamical method by about 0.01 fm~3[8}. 


The nuclear symmetry energy plays a dominant role 
in accurately describing the crust-core interface of NSs. 
It has been known that the crust-core transition density 
is highly sensitive to the isospin dependence of the nu- 
clear equation of state (EOS)[1]. Especially, the slope 
L and curvature Ksym of the nuclear symmetry energy 
as well as the L-Ksym correlation have been reported 
to be strongly correlated with the crust-core transition 
density[12, 15, 16]. In a recent work[15], a Bayesian ap- 
proach is used to infer the distribution of p; based on the 
low density constraints for neutron and symmetric nucle- 
ar matter from the effective field theory. However, the 
obtained p distribution and the correlations between p 
and the EOS parameters highly depend on the surface 
energy parameter. In the present work, we perform a 
Bayesian inference of the crust-core transition density in 
NSs based on the NS radius and neutron-skin thickness 
data. The dependence of the L-Ksym correlations pre- 
dicted in the literature on the posterior distributions of 
pı are discussed. 


The rest of paper is organized as follows. In the next 
section, we outline the theoretical framework including 
the thermodynamical method, the EOS meta-modeling 
method, the nuclear droplet model and the Bayesian in- 
ference approach. In the section II, we probe the crust- 
core transition density and its correlations with the EOS 
parameters via the NS radius and neutron-skin thickness 
data in the Bayesian framework. We also explore the 
effect of the L-Ksym correlation on the crust-core tran- 
sition density and pressure. A summary is given at the 
end. 


II. THEORETICAL FRAMEWORK 


A. Crust-core transition density and 
isospin-dependent parametric EOS for the core of 
NSs 


The crust-core transition density in the present work 
is estimated by adopting the thermodynamical approach 
under the condition that the energy per nucleon E(p, 6) 
in nuclear matter at nucleon density p and isospin asym- 
metry 6 = (Pn — Pp)/p is approximated by the isospin- 
parabolic expansion, 


E(p,5) ~ Eo(p) + Esym(p) : 0. (1) 


The crust-core transition point is obtained when the u- 
niform matter starts to be separated into a mixture con- 
taining single nucleons and clusters. The transition den- 
sity is specifically calculated by the vanishing effective 
incompressibility of uniform NS matter under 8 equilib- 
rium and charge neutrality conditions, 
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That is, one can obtain the transition density p by solv- 
ing the equation of K,, = 0. It is noteworthy to mention 
that the crust-core transition density can be overestimat- 
ed when the parabolic approximation (Eq. (1)) is used[8]. 
As indicated in Fig. 5 of Ref. [8], this overestimation 
mainly appears in the region where the parameters L 
and Ksym are large, i.e., L > 60 MeV and Kgym > -100 
MeV. However, after filtered by the NS radius data, the 
probability that they fall in the abovementioned interval- 
s is very small according to our earlier calculations[17]. 
Therefore, employing the parabolic approximation will 
hardly change the present results. Furthermore, It also 
can reduce the discrepancy between the thermodynami- 
cal method used in the present work and other approach- 
es like the dynamical method when these EOS parame- 
ters are confronted with the NS radius data. 

The transition pressure can be approximately written 
as[18] 
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where 6; is isospin asymmetry corresponding to py. Eo 
and Esym in Eqs. (1) and (2) are the energy per particle 
in symmetric nuclear matter (SNM) and the nuclear sym- 
metry energy, respectively. They can be parameterized 


as 
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where Eo(po)=-15.9 MeV is the energy per 
particle at the saturation density pọ for SN- 
M. Ko = 9%ðEo(p)/3P]- and Jo = 


2703[0° Eo(p)/3p?]|p—p represent the incompress- 
ibility and the skewness parameters of SNM at po, 
respectively. Esym(P0), L = 3p0 OE sym(P)/OP||p=p0: 
Ksym = — 9960? Esym(P)/3P’]lp= and Jsym = 
2703 [0? Esym(p)/Op*]|p=po are, respectively, the magni- 
tude, slope, curvature and skewness of nuclear symmetry 
energy at po. According to the systematics of terrestrial 
nuclear experiments and predictions of various nuclear 
theories, Ko, Esym(po), and L have the ranges of 220~ 
260 MeV, 28.5~ 34.9 MeV and 30~ 90 MeV[19-23], 
respectively. While the parameters Jo, Ksym, and 
Jsym Characterizing the nuclear EOS at high densities 
have the wide ranges of —800 < Jọ < 400 MeV, 
—400 < Ksym < 100 MeV, and —200 < Jsym < 800 MeV 
only from the theoretical predictions|24, 25]. 

In the framework of the minimum NS model, the non- 
rotating NS consists of neutrons, protons, electrons, and 
muons under the 6 equilibrium and charge neutrality 
conditions, and the relation between pressure and nu- 
cleon density in the core of NSs 
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is controlled by the energy density e(p, ô) = p[E(p,6) + 
Myn] + «1(p,6) with My and ¢;(p,6) being the average 
nucleon mass and the energy densities for leptons which 
are calculated by the noninteracting Fermi gas model|[26], 
respectively. ô can be obtained through the charge neu- 
trality condition pp = pe +p, and the $-equilibrium con- 
dition Mn — Hp = He = Hu © 40Esym(p), where u de- 
notes the chemical potential calculated by the expression 
Li = Oe(p,5)/Op; for the ith particle. 

One can construct the EOS for the core of NSs in term 
of the expressions of (4), (5) and (6). Below the crust- 
core transition density, the NV EOS(27] for the inner 
crust and the BPS EOS{1] for the outer crust are em- 
ployed. 


B. Neutron skin and nuclear droplet model 


The neutron skin thickness of a finite nucleus in the 
nuclear droplet model(DM) is obtained according to the 
expression|3, 28] 
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where e?Z/70Esym(po) is a correction term due to the 
Coulomb interaction. R = ro A!/3 is the nuclear radius. 
bn and bp denote the surface widths of the neutron and 
proton density profiles, respectively, and usually b, = 
bp = 1 fm is used in the standard DM. The quantity t in 
Eq. (7) is calculated as[3] 
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Here p4 = 0.1 fm~? for the nucleus 7°°Pb and pa = 0.08 
fm? for the nucleus 48Ca. 


C. Bayesian inference approach 


The Bayesian theorem is expressed as 
P(M|D) = CP(D|M)P(M). (10) 


Here C is a normalization constant. P(M|D) stands for 
the obtained posterior probability distribution function 
(PDF) of the model M when the data set D is given. 
P(D|M) is a likelihood function obtained through con- 
fronting the theoretical results from the model M with 
the data D, and P(M) is the prior probability of the 
model M representing the knowledge on the theoretical 
parameters of M before being compared with the data 
D. 

In the present work, we randomly sample the transi- 
tion density in the range of 0.03 fm 7? < p, < 0.2 fm ~3, 
and the matched six parameter set can be found by solv- 
ing the equation of K,,=0 by using Eq. (2). We adopt 
two ways to generate the posterior PDF of the transi- 
tion density. One of them is only based on the observed 
data of NS radii up to now, and the other is based on 
both the NS radius data and the data of the neutron- 
skin thickness. For the first way, we use the matched six 
parameters to construct the NS EOS in the framework of 
the minimum NS model as described above and put them 
in the TOV equation to compute the theoretical values 
of NS radii. Then we can obtain the likelihood of the 
transition density or the matched set of parameters by 
confronting these theoretical values of NS radii with the 
observed values using the following likelihood function, 
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g represent the observed values and its lo error bar for 
NS radii. The NS EOS generated above should satisfy 
the thermodynamical stability condition (dP/de > 0, P 
denotes the pressure inside NSs) and the causality condi- 
tion (0 < v? < c’, vs and c denote the speed of sound and 


light, respectively) at all densities, and should be stiff e- 
nough to support the maximum mass of the NS observed 
up to now. A sharp cut-off at 1.97 Mọ is used in the 
present analysis. 


For the second way, after obtaining the matched six 
parameter set using Eq. (2) through sampling the tran- 
sition density, we let them as the input into the DM to 
calculate the theoretical values of the neutron skin thick- 
ness for ?°°Pb and 48Ca. Then we can discard these 
parameter sets and the transition densities when the cal- 
culated values for the neutron skin thickness are far from 
the experimental values using the likelihood function 
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where ARD and AR}? are respectively the theoretical 
and experimental values, and o’ represent the lo error 
bar for the experimental data. After this, the rest of the 
parameter sets will be used to construct the NS EOS and 
calculate the posterior PDFs of the transition density as 
described above. 


The observed data for the NS radii and the experi- 
mental data for the neutron-skin thickness used in this 
work are summarized in Table I. These NS data in- 
clude: (i) Ria = 11.9714 km at 90% confidence lev- 
el (CFL) by analyzing the GW170817 source report- 
ed by the LIGO/Virgo Collaboration[29]; (ii) Ri4 = 
10.8776 km at 90% CFL by analyzing the same source 
of GW170817[30]; (iii) Ry.4 = 11.77}; km at 90% CFL 
reported earlier by analyzing the quiescent low-mass X- 
ray binaries (QLMXBs)[31]; (iv) Reported by the NICER 
Collaboration[33, 36] at 68% CFL, R = 12.7171 45 with 
mass of 1.34*0-13Mo([32] and R = 13.021125 with mass 
of 1.4470-13Mo[36] for the source of PSR J0030+0451, 
R = 13.7178 with mass of 2.0879-02Mo for the source of 
PSR J0740+6620. The data for the neutron-skin thick- 
ness are ARnp = 0.121+9:026 fm and ARnp = 0.283+9:971 
fm for 48Ca and 2°8Pb, respectively, which taken from 
Refs. [34, 35] recently reported by the CREX and PREX- 
2 Collaborations. 


The Metropolis-Hastings algorithm within a Markov- 
Chain Monte Carlo (MCMC) approach is adopted to gen- 
erate the posterior PDFs of the model parameters. Not 
only the PDFs of the individual parameters but the PDFs 
for the two-parameter correlations can be calculated by 
integrating over all other parameters using the marginal 
estimation approach. The initial samples in the so-called 
burn-in period are thrown away[37] so that the MCMC 
process can start from an equilibrium distribution. The 
40,000 steps for burn-in progress and the rest one million 
steps for calculating the PDF of the transition density 
are taken in the present analysis. 


TABLE I: Data for the NS radius and the neutron-skin thickness used in the present work. 


Mass(Mo) Radius R (km) Source and Reference 
1.4 11.9+1:4(90% CFL) GW170817 [29] 
1.4 10.8776 (90% CFL) GW170817 [30] 
1.4 11.7417 (90% CFL) QLMXBs [31] 
1.344915 19.71+1-14 (68% CFL) PSR J0030+0451 [32] 
1.444973 13.0475 (68% CFL) PSR J0030+0451 [33] 
2.08¢8-07  13.7+26 (68% CFL) PSR J0740+6620 [33] 
Nucleus ARnp (fm) Source and Reference 
Ca 0.12170 -628 (68% CFL) CREX([34] 

208 Ph 0.283+9:071 (68% CFL) PREX-2[35] 


III. RESULTS AND DISCUSSIONS 


A. Exploring the crust-core transition density via 
the NS observations 


The posterior PDFs of the crust-core transition density 
pı and the corresponding transition pressure P, as well as 
their correlations with the EOS coefficients are plotted in 
Fig. 1. In the calculations, the two types of prior for p 
are adopted. The first form is the uniform distribution 
which is a better choice as we have no idea about p;, and 
this case is displayed in Fig. 1 by the black curves. The 
second form is the Gaussian distribution with a average 
value of 0.078 fm~? and a standard deviation of 0.04 as 
used in Ref. [3], which is indicated by the purple curve 
in Fig. 1. In the panels located in the upper two rows we 
show the posterior PDFs of the correlations mentioned 
above using the uniform priors for p+, and these panels in 
the bottom two rows show the results using the Gaussian 
priors. These results are only based on the observed data 
of NS radii summarized in Table I. 

A two-humped posterior distribution for p is observed 
for both the uniform and Gaussian priors used in the 
calculations. The first peak is located at p, = 0.08 fm~3 
which is often used as its fiducial value in the literature, 
and the second peak is at p = 0.1 fm~3. The calculated 
68% and 90% credible intervals by using the highest pos- 
terior density interval approach for p, and P, are listed in 
Table II. We can see that, relative to the prior distribu- 
tions, the posterior PDFs of p, are narrowed down to the 
small intervals, which reflects the fact that the crust-core 
transition density is sensitive to the NS radius. It has a 
larger probability to fall in the region where p+ has values 
larger than 0.1 fm~? when a uninformative prior is used. 
This can be attributed to the correlations between p, and 
some EOS parameters and will be discussed later. 

A better constraint on p+ is found when the Gaussian 
prior is used in comparison with these using the uniform 
prior. This is easy to understand that more information 
is available for the Gaussian prior than these results from 
the uniform prior before confronting them with the NS 


radius data. The generated ranges of P,, i.e., 0.051323 


MeV /fm? using the uniform prior and 0.11928 MeV/fm? 
using the Gaussian prior at 68% confidence level as list- 
ed in Table II, can cover these calculated by using the 
meta-modelling, dynamical and thermodynamical mod- 
els, see the Table I in Ref. [15], however, the most prob- 
able values are smaller than them. Our results for p, are 
consistent with these in the literature. 


Let’s turn to explore the correlations among p+, P; and 
the EOS parameters. Here we do not concern the corre- 
lations among the EOS parameters which are consistent 
with these reported in our earlier publications[17, 40]. 
The low-order parameters such as Esym(po0) and Ko are 
hardly correlated with the transition. However, as men- 
tioned in Ref. [12], the condition K, = 0 requires a 
larger p+ as Ko increases. This phenomenon is displayed 
in Fig. 1, namely a weak positive correlation is found 
between p; and Ko for the cases using both the uniform 
and Gaussian priors. 


Consistent with the results reported in Ref. [16], the 
strong correlations among p+ and the isovector compress- 
ibility Ksym, skewness Jgym are discovered. A negatively 
correlation between p; and Ksym is inconsistent with the 
results in Ref. [16], in which the EOS parameters are 
filtered by the predictions from the effective field theo- 
ry and surface coefficients are determined by the nucle- 
ar masses in the framework of extended Thomas Fermi 
approximation method. The positively correlations be- 
tween p; and Jsym, and p+ and P, showed in Fig. 1 are in 
agreement with these in Ref. [15, 16]. The transition is 
not affected by the skewness of symmetric nuclear matter 
Jo. 

L shows a negatively correlation with p+. Interestingly, 
a positively correlation appears in the region at p, > 0.1 
fm~? for the case using the uniform prior. This does not 
happen for the case using the Gaussian prior. Do this re- 
late to the the shoulder indicated in posterior PDF of p 
using the uniform prior in Fig. 1 ? To answer this ques- 
tion, we plot the L-Kgym correlations from three types of 
calculations as indicated in Fig. 2, i.e., using the uniform 
and Gaussian priors based on the NS radius data, using 
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FIG. 1: (color online)Posterior probability distribution functions of the crust-core transition density and pressure as well as 
their correlations with the EOS parameters. The calculations are performed by adopting the uniform (shown by blue curves) 
and Gaussian (shown by red curves) priors for the transition density based on the NS radius data. For comparison, the uniform 
(shown by black curve) and Gaussian (shown by purple curve) forms for the prior distributions of the transition density are 


included. The correlations in the upper (bottom) two rows are calculated by using the uniform (Gaussian) prior for the 
transition density. 


TABLE II: Most probable values and the 68% and 90% credible intervals of p: and P, calculated by using the uniform prior, 
the Gaussian prior, the neutron-star plus neutron-skin thickness data (NS+NST), the L— Ksym correlation in Ref. [38] by Holt 
et al., the L — Ksym correlation in Ref. [39] by Mondal et al., the L — Ksym correlation in Ref. [24] by Tews et al.. 


Uniform prior Gaussian prior NS+NST Holt et al.(2018) Mondal et al.(2017) Tews et al. (2017) 
pi 0.0875'995:0-085 5:25 0.0875'93,0-084 5:93” 0.07545 '91 50-075 25:91, 0-085 5'095,0-087 5.015 0-1t001 0.14002 0.155:097,0-17 5.02 
P, We p05 so Ud a 0.179758,0-1 75:68  0.164006:0-16Ż008 0.6Ż0:14,0.6t40.28 —_(0.61514,0.6 15°58 
the uniform prior based on both the NS radius and neu- correlated to the EOS parameters. 


tron skin thickness data. There are two phenomena in 
Fig. 2 we can see. One is the anti-correlations shown in 
left and right panels, and highly weak correlation shown 
in the middle panel between L and Ksym, and the other 
is, as shown in the left panel in Fig. 2, Ksym has a high 
probability to stay in the region where Ksym is extreme- 
ly negative, i.e. Ksym < -200 MeV. The latter is obvi- 
ously responsible for the shoulder because the shoulder 
for the posterior PDF of p disappears as shown in Fig. 
3 though the L-Ksym anti-correlations are same in their 
calculations in the left and right panels. Therefore, Ksym 
plays a more important role in constraining the crust-core 
transition density of NSs than the L-Ksym correlation as 
studied in Ref. [12]. The transition pressure are weakly 


B. Effect of neutron-skin thickness and comparison 
with other calculations 


It has been known that the neutron-skin thickness is a 
effective probe to the nuclear symmetry energy, especially 
its slope parameter[41, 42]. The latter plays an impor- 
tant role in determining the crust-core transition density 
of NSs. Shown in Fig. 3 are the results same as these 
in Fig. 1 but based on not only the NS radius data but 
also the neutron-skin thickness data. In the calculations, 
we first perform a Bayesian inference of the coefficients 
Esym(po), L and Ksym in the framework of the nucle- 
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FIG. 2: (color online)Correlations between L and Kym calculated from the NS radius data by using the uniform (Left) and 
Gaussian (Middle) priors, and from both the neutron-star radius and neutron-skin thickness data (Right) by using the uniform 
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FIG. 3: (color online)Posterior probability distribution functions of the crust-core transition density and pressure as well as 


their correlations with the EOS parameters. 


In the calculations, both the neutron-skin thickness and NS radius data are 


adopted. The prior distribution (shown by black curve) for the transition density is included. 


ar droplet model as described in Sec. IIB based on the 
neutron-skin thickness data reported by the CREX[34] 
and PREX-2[35] Collaborations as listed in Table I. Af- 
ter that, we regard the obtained distributions of these 
parameters and the matched transition density as their 
priors to infer the posterior PDF of the transition den- 
sity based on the NS radius data within the minimum 
NS model. It is worth noting that there exists significant 
controversy about constraining the slope of the symmetry 
energy through the neutron-skin thickness data of *%Ca 
and ?°8Pb reported by the CREX and PREX-2 Collabo- 
rations. A recent study demonstrates that the ranges for 
the slope L obtained from the CREX and the PREX-2 
are completely inconsistent[43]. More accurate measure- 
ments are needed. Fortunately, there are many methods 
to determine the neutron-skin thickness, i. e., the con- 
figurational information entropy analysis[44, 45]. 


As seen in Fig. 3, the correlations are roughly same as 


these in Fig. 1. A weak anti-correlation between the sym- 
metry energy magnitude Esym(po) and pz is observed due 
to the constraint from the neutron-skin thickness data on 
Esym(po). The range for Ksym is smaller than these us- 
ing only the NS radius data, which is also because of the 
effect of the neutron-skin thickness data. The constrain- 
t on pz obtained is improved in comparison with these 
only based on the NS data as listed in Table II. This 
is because the parameters L and Ksym are better con- 
strained when the neutron-skin thickness data are taken 
into consideration. 


In Fig. 4 we compare our results with these inferred us- 
ing a compressible liquid-drop model within a Bayesian 
framework[15], in which there are two filters, i.e., the 
low density (LD) behavior of the energy functionals that 
should be rigorously limited in the uncertainty interval- 
s of the effective field theory calculation for symmetric 
and pure neutron matter in Ref. [46] and the high den- 
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FIG. 4: (color online)Posterior probability distribution func- 
tions of the crust-core transition density and pressure based 
on only the NS radius data by using the uniform (shown by 
blue curve) and Gaussian (shown by red curve) priors, and 
based on both the NS radius and neutron-skin thickness data 
by using the uniform prior indicated by NS+NST (shown by 
black curve). For comparison, the results in Ref. [15] shown 
by LD-EPJA2019 and HD-EPJA2019 are included. 


sity (HD) behavior of the functionals that should obey 
the conditions such as positive symmetry energy at all 
densities, causality and so on[15]. One can see that our 
results overlap with theirs. There is a long tail for P, 
for the case using the uniform prior in the calculations, 
which is because the large probability of p, exists in the 
region at p; > 0.1 fm~3. Our most probable values for 
P, are smaller than these from Ref. [15]. 


C. Effect of L-Ksym correlations 


It has been reported that L-K sym correlations have a 
significant impact on both the crust-core transition den- 
sity and pressure[12]. To further explore this effect, we 
take three typical L-Ksym correlations predicted in the 
literature as the priors to infer the posterior PDF of p+. 
For completeness and ease of discussions, here we de- 
scribe briefly the three correlations. The first one is based 
on the theoretical predictions from 240 Skyrme Hartree- 
Fock and 263 relativistic mean-field calculations[39] by 


„ MeV] 


sy 


K 


— Mondal et al. (2017) 
— Holt et al. (2018) 
— Tews et al. (2017) 


40 45 50 55 60 65 70 75 30 
L [MeV] 


FIG. 5: (Color online) L — Ksym correlations from Tews et 
al. [24] (shown by red curves), Mondal et al. [39] (shown by 
black curves) and Holt et al. [38] (shown by green curves). 
The solid and dashed lines respectively denote the means and 
the boundaries of each correlation. Taken from Ref. [12]. 


Mondal et al. and written as 


Kym = (—4.97£0.07)(3Esym(po)—L)-+66.80+2.14 MeV. 

(13) 
The second one is taken from Ref. [24] by Tews et al. 
and written as 


Ksym = 3.50L — 305.67 + 24.26 MeV. (14) 


In the framework of the Fermi liquid theory, Holt and 
Lim[38] derived the expressions L = 6.70Esym(po) — 
148.60 +4.37 MeV and Ksym = 18.50Esym(po) — 613.184 
9.62 MeV. Thus one can formulate the last relation be- 
tween L and Ksym as[12] 


Ksym = 2.76L — 203.07 + 21.69 MeV. (15) 


These correlations are shown in Fig. 5. The discrepancy 
among the correlations by Mondal et al. and Tews et al. 
is not too much larger because they are from the same 
sets of theory predictions. However, These depicted in 
Figs. 2 and 5 are highly different. This illustrates that 
the L-Ksym correlations are strongly model dependent 
and constraining them has a long way to go. 

In the Bayesian inference approach, after considering 
the abovementioned correlations, L and Ksym will not 
be independent when we randomly sample them between 
their specific ranges. In the calculations performed in this 
subsection, the uniform prior for p; and only the NS ra- 
dius data are employed. The generated posterior PDFs 
for p, and P, are represented in Fig. 6, and the corre- 
sponding confidence intervals are summarized in Table 
II. One can see that, as stated in Ref. [12], the L— Ksym 
correlations play a significant role in constraining both 
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FIG. 6: (color online)Posterior probability distribution functions of the crust-core transition density and pressure as the three 
correlations of L — Ksym as in Fig. 5 are used in the calculations. 


the transition density and pressure. There are at least 
two observations we can obtain: (i)The results from the 
relations by Mondal et al. and Tews et al. are completely 
consistent. This is mainly because this two correlation- 
s largely overlap within the allowed error limits; (ii)The 
most probable values from the relations by Holt et al. are 
highly different from the other two cases. For the latter, 
the transition density and pressure have the larger values 
than these for the former. This is due to the fact that 
Ksym is not extremely negative in the relation by Holt et 
al., i.e. Ksym having the values higher than about -115 
MeV as shown by green curves in Fig. 5. The present 
results are consistent with these in Ref. [12], in which the 
authors studied the effects of L — Ksym correlations on 
the crust-core transition density and pressure by adopt- 
ing the fixed values of the other parameters in Eqs. (4) 
and (5). 


IV. SUMMARY 


In summary, using the thermodynamical approach to 
calculate the crust-core transition density and an explic- 
itly isospin-dependent parametric EOS for the core of 
NSs within the minimum NS model, the Bayesian in- 
ference of the crust-core transition density based on the 
NS radius and neutron-skin thickness data is performed. 
The uniform and Gaussian forms for prior distributions of 


[1] G. Baym, C. Pethick, P. Sutherland, The ground s- 
tate of matter at high densities: equation of state 
and stellar models. The Astrophys. J. 170, 299 (1971). 
doi:10.1086/151216 

[2] N. Paar, C.C. Moustakidis, T. Marketin, et al., 
Neutron star structure and collective excitations of 


transition density are employed in the calculations. It is 
found that the transition density has a larger probability 
to take the values larger than 0.1 fm~? when the unifor- 
m prior is used, which does not happen as the Gaussian 
prior is used. To take the values smaller than -200 MeV 
for Ksym is responsible for this phenomenon. 


The negatively (positively) correlations between p; and 
L and between p; and Ksym (between p; and Ko, between 
pz and Jsym and between p; and P;) are observed. Based 
on the NS radius data reported so far, the generated 
68% confidence intervals for p, are respectively O08 poe 
fm~3 and 0.0810'}3 fm— when the uniform and Gaussian 
priors are adopted. When the EOS parameters are first 
filtered by the neutron-skin thickness data, the obtained 
one is 0.075790)" fm-°. We also check the impact of 
L-Ksym correlations on the posterior PDFs of p; and Py. 
The results from the L-Ksym relations reported by Tews 
et al. [24] and Mondal et al.[39] are completely over- 
lapped, and the results reported from Holt et al. [38] are 
smaller than these from the other two L-K sym relations. 


Acknowledgement: This work is supported by the 
Shanxi Provincial Foundation for Returned Overseas 
Scholars under Grant No 20220037, the Natural Sci- 
ence Foundation of Shanxi Province under Grant No 
20210302123085, and the discipline construction project 
of Yuncheng university. 


finite nuclei. Phys. Rev. C 90, 
doi:10.1103/PhysRevC.90.011304 

[3] M. Centelles, X. Roca-Maza, X. Vinas, et al., Nu- 
clear symmetry energy probed by neutron skin thick- 
ness of nuclei. Phys. review letters 102, 122502 (2009). 
doi:10.1103/PhysRevLett.102.122502 


011304 (2014). 


[4] 


[5] 


[6] 


[8] 


(9 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


C. Horowitz, J. Piekarewicz, Neutron star structure and 
the neutron radius of p 208 b. Phys. Rev. Lett. 86, 5647 
(2001). doi:10.1103/PhysRevLett.86.5647 

L. Tsaloukidis, C. Margaritis, C.C. Moustakidis, Effects 
of the equation of state on the core-crust interface of 
slowly rotating neutron stars. Phys. Rev. C 99, 015803 
(2019). doi:10.1103/PhysRevC.99.015803 

J. Fang, H. Pais, S. Pratapsi, et al., Crust-core transition 
of a neutron star: Effects of the symmetry energy and 
temperature under strong magnetic fields. Phys. Rev. C 
95, 062801 (2017). doi:10.1103/PhysRevC.95.062801 

C. Gonzalez-Boquera, M. Centelles, X. Viñas, et al., 
Core-crust transition in neutron stars with finite-range 
interactions: The dynamical method. Phys. Rev. C 100, 
015806 (2019). doi:10.1103/PhysRevC.100.015806 

J. Xu, L.W. Chen, B.A. Li, et al., Nuclear constraints on 
properties of neutron star crusts. The Astrophys. J. 697, 
1549 (2009). doi:10.1088/0004-637X /697/2/1549 

J.M. Lattimer, Y. Lim, Constraining the symmetry pa- 
rameters of the nuclear interaction. The Astrophys. J. 
771, 51 (2013). doi:10.1088/0004-637X/771/1/51 

S.S. Bao, H. Shen, Impact of the symmetry ener- 
gy on nuclear pasta phases and crust-core transition 
in neutron stars. Phys. Rev. C 91, 015807 (2015). 
doi:10.1103/PhysRevC.91.015807 

C.C. Moustakidis, Effect of the symmetry ener- 
gy on the location of the inner edge of the neu- 
tron star crust. Phys. Rev. C 86, 015801 (2012). 
doi:10.1103/PhysRevC.86.015801 

B.A. Li, M. Magno, Curvature-slope correlation of nu- 
clear symmetry energy and its imprints on the crust- 
core transition, radius, and tidal deformability of canon- 
ical neutron stars. Phys. Rev. C 102, 045807 (2020). 
doi:10.1103/PhysRevC.102.045807 

Z.W. Liu, Z. Qian, R.Y. Xing, et al., Nuclear fourth-order 
symmetry energy and its effects on neutron star proper- 
ties in the relativistic hartree-fock theory. Phys. Rev. C 
97, 025801 (2018). doi:10.1103/PhysRevC.97.025801 

H. Pais, A. Sulaksono, B.K. Agrawal, et al., Corre- 
lation of the neutron star crust-core properties with 
the slope of the symmetry energy and the lead 
skin thickness. Phys. Rev. C 93, 045802 (2016). 
doi:10.1103/PhysRevC.93.045802 

T. Carreau, F. Gulminelli, J. Margueron, Bayesian anal- 
ysis of the crust-core transition with a compressible 
liquid-drop model. The Eur. Phys. J. A 55, 188 (2019). 
doi:https://doi-org/10.1140/epja/i2019-12884-1 

S. Antić, D. Chatterjee, T. Carreau, et al., Quantifying 
the uncertainties on spinodal instability for stellar matter 
through meta-modeling. J. Phys. G: Nucl. Part. Phys. 46, 
065109 (2019). doi:10.1088/1361-6471/abla51 

W.J. Xie, B.A. Li, Bayesian inference of high-density 
nuclear symmetry energy from radii of canonical 
neutron stars. The Astrophys. J. 883, 174 (2019). 
doi:10.3847 /1538-4357 /ab3f37 

J.M. Lattimer, M. Prakash, Neutron star obser- 
vations: Prognosis for equation of state con- 
straints. Phys. reports 442, 109-165 (2007). 
doi:https://doi.org/10.1016/j.physrep.2007.02.003 

P. Danielewicz, R. Lacey, W.G. Lynch, Determination of 
the equation of state of dense matter. Sci. 298, 1592- 
1596 (2002). doi:10.1126/science.1078070 

S. Shlomo, V. Kolomietz, G. Colo, Deducing the 
nuclear-matter incompressibility coefficient from 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


[32] 


[33] 


[34] 


[35] 


The Eu- 
(2006). 


data on isoscalar compression modes. 
r. Phys. J. A-Hadrons Nucl. 30, 23-30 
doi:https://doi.org/10.1140/epja/i2006-10100-3 
J. Piekarewicz, Do we understand the incompressibility 
of neutron-rich matter? J. Phys. G: Nucl. Part. Phys. 
37, 064038 (2010). doi:10.1088/0954-3899/37/6/064038 
B.A. Li, X. Han, Constraining the neutron—proton ef- 
fective mass splitting using empirical constraints on the 
density dependence of nuclear symmetry energy around 
normal density. Phys. Lett. B 727, 276-281 (2013). 
doi:https://doi.org/10.1016/j.physletb.2013.10.006 

M. Oertel, M. Hempel, T. Klähn, et al., Equations of 
state for supernovae and compact stars. Rev. Mod. Phys. 
89, 015007 (2017). doi:10.1103/RevModPhys.89.015007 
I. Tews, J.M. Lattimer, A. Ohnishi, et al., Symmetry 
parameter constraints from a lower bound on neutron- 
matter energy. The Astrophys. J. 848, 105 (2017). 
doi:10.3847 /1538-4357 /aa8db9 

N.B. Zhang, B.J. Cai, B.A. Li, et al., How tightly 
is the nuclear symmetry energy constrained by a u- 
nitary fermi gas? Nucl. Sci. Tech. 28, 1-7 (2017). 
doi:https: //doi.org/10.1007/s41365-017-0336-2 

N.B. Zhang, B.A. Li, J. Xu, Combined constraints 
on the equation of state of dense neutron-rich matter 
from terrestrial nuclear experiments and observation- 
s of neutron stars. The Astrophys. J. 859, 90 (2018). 
doi:10.3847 /1538-4357 /aac027 

J.W. Negele, D. Vautherin, Neutron star matter at sub- 
nuclear densities. Nucl. Phys. A 207, 298-320 (1973). 
doi:https: //doi.org/10.1016/0375-9474(73)90349-7 

M. Warda, X. Vinas, X. Roca-Maza, et al., Neu- 
tron skin thickness in the droplet model with surface 
width dependence: Indications of softness of the nucle- 
ar symmetry energy. Phys. Rev. C 80, 024316 (2009). 
doi:10.1103/PhysRevC.80.024316 

B.P. Abbott, R. Abbott, T. Abbott, et al, G- 
w170817: Measurements of neutron star radii and equa- 
tion of state. Phys. review letters 121, 161101 (2018). 
doi:10.1103/PhysRevLett.121.161101 

S. De, D. Finstad, J.M. Lattimer, et al., Tidal deforma- 
bilities and radii of neutron stars from the observation 
of gw170817. Phys. review letters 121, 091102 (2018). 
doi:10.1103/PhysRevLett.121.091102 

J.M. Lattimer, A.W. Steiner, Constraints on the sym- 
metry energy using the mass-radius relation of neu- 
tron stars. The Eur. Phys. J. A 50, 1-24 (2014). 
doi:https://doi.org/10.1140/epja/i2014-14040-y 

T.E. Riley, A.L. Watts, S. Bogdanov, et al., A nicer 
view of psr j0030+ 0451: millisecond pulsar parameter 
estimation. The Astrophys. J. Lett. 887, L21 (2019). 
doi:10.3847 /2041-8213/ab481c 

E. Fonseca, H.T. Cromartie, T.T. Pennucci, et al., Re- 
fined mass and geometric measurements of the high- 
mass psr j0740+ 6620. The Astrophys. J. Lett. 915, L12 
(2021). doi:10.3847 /2041-8213 /ac03b8 

D. Adhikari, H. Albataineh, D. Androic, et al., Pre- 
cision determination of the neutral weak form fac- 
tor of ca 48. Phys. Rev. Lett. 129, 042501 (2022). 
doi:10.1103/PhysRevLett.129.042501 

D. Adhikari, H. Albataineh, D. Androic, et al., Ac- 
curate determination of the neutron skin thickness 
of pb 208 through parity-violation in electron s- 
cattering. Phys. review letters 126, 172502 (2021). 
doi:10.1103/PhysRevLett.126.172502 


[36] 


[37] 


[38] 


[39] 


[40] 


[41] 


[42] 


[43] 


[44] 


[45] 


[46] 


M.C. Miller, F. Lamb, A. Dittmann, et al., The 
radius of psr j0740+ 6620 from nicer and xmm- 
newton data. The Astrophys. J. Lett. 918, L28 (2021). 
doi:10.3847/2041-8213/ac089b 

R. Trotta, Bayesian methods in cosmol- 
ogy. arXiv preprint arXiv:1701.01467 
doi:https: //doi.org/10.48550/arXiv.1701.01467 

J.W. Holt, Y. Lim, Universal correlations in 
the nuclear symmetry energy, slope parameter, 
and curvature. Phys. Lett. B 784, 77-81 (2018). 
doi:https://doi.org/10.1016/j.physletb.2018.07.038 

C. Mondal, B. Agrawal, J. De, et al., Interdependence 
of different symmetry energy elements. Phys. Rev. C 96, 
021302 (2017). 

W.J. Xie, B.A. Li, Bayesian inference of the symmetry 
energy of superdense neutron-rich matter from future ra- 
dius measurements of massive neutron stars. The Astro- 
phys. J. 899, 4 (2020). doi:10.3847/1538-4357/aba271 
J. Xu, Constraining isovector nuclear interactions with 
giant dipole resonance and neutron skin in 208pb from a 
bayesian approach. Chin. Phys. Lett. 38, 042101 (2021). 
doi:10.1088/0256-307X/38/4/042101 

J. Xu, W.J. Xie, B.A. Li, Bayesian inference of nuclear 
symmetry energy from measured and imagined neutron 
skin thickness in sn 116, 118, 120, 122, 124, 130, 132, 
pb 208, and ca 48. Phys. Rev. C 102, 044316 (2020). 
doi:10.1103/PhysRevC.102.044316 

S. Tagami, T. Wakasa, M. Yahiro, Slope parameters de- 
termined from crex and prex2. Results Phys. 43, 106037 
(2022). doi:https://doi.org/10.1016/j.rinp.2022.106037 
H.L. Wei, X. Zhu, C. Yuan, Configurational infor- 
mation entropy analysis of fragment mass cross dis- 
tributions to determine the neutron skin thickness 
of projectile nuclei. Nucl. Sci. Tech. 33, 111 (2022). 
doi:https: //doi.org/10.1007/s41365-022-01096-w 

C.W. Ma, Y.P. Liu, H.L. Wei, et al., Determina- 
tion of neutron-skin thickness using configurational 
information entropy. Nucl. Sci. Tech. 33, 6 (2022). 
doi:https://doi.org/10.1007/s41365-022-00997-0 

C. Drischler, K. Hebeler, A. Schwenk, Asymmet- 
ric nuclear matter based on chiral two- and three- 
nucleon interactions. Phys. Rev. C 93, 054314 (2016). 
doi:10.1103/PhysRevC.93.054314 


10 


